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Abstract 

We investigate the dynamics of entanglement and quantum discord of two qubits in liquid state homonuclear NMR. 
Applying a phenomenological description for NMR under relaxation process, and taking a group of typical parameters 
of NMR, we show that when a zero initial state |00) experiences a relaxation process, its entanglement disappears 
completely after a sequence of so-called sudden deaths and revivals, while the quantum discord retains remarkable 
values after a sequence of oscillations. That is to say, the quantum discord is more robust than entanglement. 
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1. Introduction 



Entanglement is one of the most striking features for our deeper understanding of the quantum world, and is widely 
seen as the main resource for quantum information processing [ However, entanglement is not the only type of 
useful correlation. It has been found that entanglement is not necessary for deterministic quantum computation with 
one pure qubit (DQC1) [3—5], but another quantum correlation, introduced by |@,|3], called quantum discord (QD), 
is responsible for the computational efficiency of DQC1. Besides, QD has also been used for some studies such as 
quantum phase transition and the process of Grover searchfl. Recently, the dynamics of QD has attracted more 
■ attentions, such as the robustness of QD to sudden death [)To|], QD in Heisenberg models | IT] and Non-Markovian 
effect of QD [HI. 

Nuclear magnetic resonance (NMR) has been used to demonstrate basic concepts of quantum information pro- 
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\Q , cessing, and has been leading the field both in terms of number of qubits and control techniques 013111411 . The studies 
of dynamics characteristics and the implementations for quantum information processing of NMR systems have also 



received a great deal of attentions in recent years [ 15-20 1. 



Due to the fundamental and practical significance, we wish to investigate the dynamics behaviors of entanglement 
and QD of two qubits in NMR systems. The rest of this paper is organized as follows. In Sec. 2, we present a brief 
review of the entanglement, QD and liquid state two qubits NMR systems. In Sec. 3, we investigate the dynamics of 
entanglement and quantum discord of two qubits NMR without relaxation. In Sec.4, applying a phenomenological 
method, we investigate the dynamics of entanglement and quantum discord of two qubits NMR with relaxation. Sec. 5 
is a brief summary. 

2. Entanglement,quantum discord, NMR systems 

We first recall some basics about entanglement, QD and NMR. 
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2.1. Concurrence of two qubits states 

Two quantum systems A and B, described by the Hilbert spaces H A and H B , and the composite system AB is 
described by the Hilbert space H A <g> H B . A state p on H A ® H B is called a separable state if and only if it can be 
written as the form 



A ^ „B 



P?, (1) 



where n is a positive integer, > 0, Yli=\ Pi - 1> \pf }" = i are density operators on H A , {pf }" =1 are density operators 
on H B . A state is called an entangled state if it is not a separable state. For two qubits, one of the most widely used 
entanglement measures is the so-called concurrence, we will use it as the quantifier in this paper. For two qubits state 
p, the concurrence of p is defined as 112 111 

C(p) = max{0, h-h.-h~ M), (2) 

where A\, A2, /I3, A4 are the square roots of the eigenvalues of the matrix pp in decreasing order, and 

p = (o-y <g> o- v )p*(cr v <g> o- v ), (3) 

cr y is the Pauli matrix, p* is the complex conjugate of p. For convenience, in Eq.(2) we denote 

C(p) = Ai-A 2 -A 3 - A 4 , (4) 

and call it pseudo-concurrence. 

2.2. Quantum discord of two qubits states 

The quantum discord of p on H A ® H B (with respect to system A) is defined as [6, 7] 

D A (p) = S(p A )-S(p AB )+ inf Yp a S(p B ). (5) 

In Eq.(5), S(-) is the entropy function, that is, e.g., ^(p 4 ) = -fr^p 4 \og 2 p A ). n& = dimH A , inf takes all projective 
measurements {|ff)(o'|(g)/}'^ I on system A, that is, {|a)}^| =1 is an arbitrary orthonormal basis for H A . We use / to denote 
the identity operators on H A or on H B . p a = tr^ {\a)(a\ <8> Ip\a){a\ ® I), p B = \a)(a\ ® Ip\a){a\ ® // p a . 

Quantum discord captures more correlation than entanglement, this can been seen by Eq.(l) and the fact that |01 

£> A (p) = 0^p = J]^l a >< ff l®^' (6) 

a=l 

where, {Iff)}"'^ is an arbitrary orthonormal basis for H A , {p B }"* =l are density operators on H B , p a > 0, p a = 1. 
A state p satisfying Da(p) — is called zero quantum discord state or classical state. Comparing Eq.(l) with Eq.(6), 
it is obvious that a classical state must be a separable state. 

Eq.(5) is difficult to optimize, even for two qubits states, up to now only few special states were found to allow 



analytical solutions 1221 12311 . A technical definition of quantum discord, called geometrical measure, was introduced 



asj2j, 

D G A {p) = inf tr[(p-cr) 2 ], (7) 

(T 

where inf takes all cr that Da(o~) = 0. Evidently, 

D A (p) = <=> D G A {p) = 0. (8) 
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One of the most elegant results about the definition of D^(p) is that it allows analytical expressions for all two qubits 
states. More specifically, for any two qubits state p which can be written as 



3 3 3 



p = hi ® I + ^ XiCTi ® / + ^ 37/ ® o-j + ^ Tijo-i ® cr ; ), (9) 



i=i j=\ ij=i 



then |2J|, 

1 3 3 

D A(P) = 4<Z x " + Z ^' " Amax) ' (10) 
i=i U=i 

Where, o"i = cr x , cr 2 = cr v , cr 3 = o~ z , are Pauli matrices, fx,}? =1 , {yj} 3 =1 , {7y}? . =1 are all real number sets. A m - dx is the 
largest eigenvalue of the matrix xx' + TT', x = (x\,X2, xj,)', T is the matrix (7^), t means transpose. Also 

Xj = tr\pcTi ® r\,yj = tr[pl ® cry], Ty = fr[pcr ; ® crj], (11) 

(The explicit expressions for {x,}^ and {Ty}? - =1 by the elements of p, see Appendix A.) 

2.3. Two qubits in liquid state NMR 

The Hamiltonian of two qubits in NMR (we focus on the liquid state NMR) is well described by [13] 

H(t) = H sys + H rf (t), (12) 
ha>\ Ka>2 

H svs - — — o- z ® / — / ® cr z + hJo- z ® cr z , (13) 

H rf (t) = — ^- [cos(wf)o" x ® / - sin( ut)cr x ® /] 2_ [ C os(tot)I ® cr x - sin(wf)/ ® cr v ] . (14) 

In H svs , oj\ = Boji, ol>2 = ^072, are called Lamor frequencies of two qubits (two nuclei), containing the chemical 
shifts. Bo is the static magnetic field along z direction, j\, j2 are gyromagnetic ratios of two qubits. Typical values 
of at\, u)\ are a few 10 8 Hz, and chemical shifts of a few 10 3 Hz to a few 10 4 Hz. J describes the spin-spin couplings 
including direct dipole-dipole coupling and indirect through-bond coupling. Spin-spin couplings are very small in 
NMR compared with the Lamor frequencies, for instance, a few 10 2 Hz. 

In H,-f(t), gi = B\yu g2 = B\j2, B\ is the applied magnetic field rotating in x-y plane at frequency u>, at or near 
the Lamor frequencies a>i, cl>2- Typical values of g\, g2 are up to 10 5 Hz. 



3. Dynamics of entanglement and quantum discord of two qubits in NMR without relaxation 

We study the dynamics of two qubits NMR without relaxation. When the evolution time is much shorter than the 
relaxation time scales, the system can be approximated as an isolated system, and it evolves obeying the Schrodinger 
equation 

m-mt)) = H(t)mt)). d5) 

Notice that H r j{t) in Eq.(14) is time dependent, in order to cancel the time t in Hamiltonian H(t), we put 

mt)} = exp[^(cr, ® / + / ® o-Mfit))- (16) 

This transformation is widely used in NMR theory, sometimes called the method of rotating frame. With this trans- 
formation, we get 
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m-w)) = my®), 

H = cr 7 



, h((02 ~ OJ) T 3- t %1 , %2 , 

i / / ® <r z + hJcr z ® cr z — <r x ® / — / ® cr*, 



mo» = i^(o)). 

In above reductions, we have used the facts that for any real number A, 

exp(iAo-j) = cos(/l)/ + i sin(/i)cr ; , 
exp(-iAo~ z )o~ x exp(iAo~ z ) — cos(2/l)cr x + sin(2A)cr y , 
expi-iAcr^o-yexpiiAcr^ = cos(2A)<r y - sm(2A)cr x . 

Since // is time independent, together with Eq.(16) and Eq.(19), we get 

\ip{i)) = exp(^i/)|^(0)>, 

W)) = WMO)), 

itot t — 

U(t) = exp[— (o% ® / + / ® <r z )] ex P(^^)- 

In matrix notation, 



exp[— (<r z ®/ + /®<r z )] 
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H = -h 



Suppose the initial state is |00), that is 
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then 



After time t, p will evolve to p(t) as 
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C(p) = 0, D c A (p) = 0. 



p(t) = U(t)pU\t), 



(28) 



(29) 



(30) 



where | means Hermitian conjugate. Then C(p(t)), D^(p(t)) can be calculated by Eq.(2) and Eq.(10). 

Fig.l shows C(p(f)) and D°(p(t)) at t = 10~ 5 s, t = IQ^s, t = 10 _1 i. Where we take the parameters of NMR as 
^p- = 3 x 10 8 ffz, i ^ = \0 4 Hz, J = 3 x 10 2 //z, f = 5 x \0 4 Hz, f = 5 X 10 4 Hz. From Fig.l we see that, only 
near the frequency "'+" 2 = 3 x 10 8 //z, entanglement and quantum discord are remarkably generated and evolved. It 
is a common phenomenon in NMR that the applied magnetic B\ can effectively influence or control the nuclear spins 
at or near the resonance frequencies. 

Fig.2 shows C(p(f)) and D°(p(t)) at the frequency = 3 x 10 8 #z, and we take the parameters of NMR as in 

Fig.l. From Fig.l and Fig.2 we see that the dynamical behaviors of entanglement and quantum discord are similar in 
the evolution where we ignore the relaxation. 
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Figure 1: Pseudo-concurrence and quantum discord at time t = 10 s, t = 10 s, t = 10 s. 
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Figure 2: Pseudo-concurrence and quantum discord at frequency a> - 
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4. Dynamics of entanglement and quantum discord of two qubits in NMR with relaxation. 

When the evolution time gets comparable to the relaxation time scales, the relaxation effects must be taken into 
account. The relaxation process in NMR can be described in a phenomenological way ( I25ll ; lfl4ll . 2.10) 



at in 



T i} = 6 {J Ti + (l-6ij)T 2 . 



(31) 

(32) 
(33) 



Where T\ is the energy relaxation rate, T2 is the phase randomization rate. Theoretical calculations and experimental 
measurements for T\ and T2 are well developed. The typical value of T\ is tens of seconds, and T 2 is easily on the 
order of 1 second or more. 
Similar to Eq.(16), we let 



iu>t iu>t 
p(t) = exp[— (cr- 0/4-7® o%)]cr(f) exp[ — ^-(°"2 ®I + I® cr z )], 



(34) 



notice that 



P(0) = cr(0). 
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(35) 



Using Eqs. (20-22), after some straightforward calculations, then Eq.(31) becomes 



— — = —[H,cr(t)] - 
at in 
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(36) 



Eq.(36) is a first-order linear nonhomogeneous ordinary differential equations with constant coefficients in the func- 
tions <T\i(t), (Tn(i), O"i3(0, o~\4(i), o"2i(f), o"44(f). If we denote <x(f) and the last matrix in Eq.(36) by 



-&{t) = (cr n (f), (Ti2(fX 0-13(0, cr I4 (f),0-2i(f),-,cr 4 4(f))', 

~?f t \ 0-12(0) -to, cri 3 (0) _ iu) , 0-14(0) _ 2 to( O~2l(0) to, 0-44(0) ,, 

/(?) = ( , e , e , e , e ), 

J T x T 2 T 2 T 2 T 2 T x 

then, Eq.(36) can be rewritten as a compact form 

^&(t)=A-&(t)+^(t\ 
at 



(37) 
(38) 



(39) 



where, A is a 16 x 16 matrix which is independent of t (the explicit expression of A, see Appendix B). So the general 
solution of Eq.(39) is 



o*(0 = exp(fA)o*(0) + f exp[(f - s)A]J(t)ds. 
Jo 



We now take the zero initial state |00) for p, hence 

o*(0) = (1,0,0,0,0,...,0)'. 

?(0 = (^,o,o,o,o,...,oy. 

1 1 



(40) 

(41) 
(42) 



Then calculate Eqs. (40,34,2, 10), we will get C(p(t)) and D^(p(t)). From the discussions about Fig.l and Fig. 2, in 
Fig. 3 below we only discuss the dynamics of C(p(t)) and D^(p(t)) at frequency 0Jl *'° 2 = 3 x 10 & Hz, and we take the 
parameters of NMR as in Fig. 1 . 

From Fig. 3 we see that, when the evolution time is much shorter than the relaxation scales, the dynamics of 
entanglement and quantum discord are similar. But for a longer time, namely, near relaxation scales, the entanglement 
manifests a sequence of sudden deaths and revivals, and finally disappears completely. At the same time, the quantum 
discord, after a sequence of oscillations, still retains remarkable values. 

In a dynamical process, the phenomenon of entanglement disappears in a finite time is named entanglement sudden 
death [26]. This is an interesting phenomenon in contrast to the usual intuition that decoherence will be in an infinite 
time. Also, this is an important phenomenon for quantum information processing since much of quantum information 
processing rely on entanglement Ill|,|2j]. 



5. Summary 

In summary, we investigated the dynamics of entanglement and quantum discord of two qubits in liquid state 
homonuclear NMR systems. We showed that the dynamical behaviors of entanglement and quantum discord are 
similar in short time where relaxation effects can be neglected. When the time is long enough to be comparable to the 
relaxation rates, the entanglement manifests the phenomenon of sudden deaths and survivals, and at last disappears 
completely. Meanwhile the quantum discord retains remarkable values. That is to say, quantum discord is more 
robust than entanglement against the relaxation processes. Hence quantum algorithms based only on quantum discord 
correlation may be more robust than those based on entanglement. 
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Appendix A. Explicit expressions for {je,-} 3 and {Ty} 3 by the elements of p in Eq.(ll) 

I — L l,J—L 

X\ = Pl3 +P24 + P31 + P42 
X2 = i(Pl3 +P24 -P31 ~P4l) 
X-i = pn +P22 -P33 -P44 
T\ 1 = Pl4 + P23 + P32 + P41 
T\2 - i(pu -P23 +P32 -P4l) 
7^13 = Pl3 - P24 + P31 - P42 
7^21 = KPU +P23 -P32 -P4l) 
7^22 = -Pl4 + P23 + P32 - P41 
7^23 = i(Pl3 -P24 -P31 +P42) 
T31 = pn + P21 ~ P34 - P43 
7^32 = i(Pl2 - P21 - P34 + P43) 
7^33 =Pll -P22 -P33 +P44 
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Appendix B. Explicit expression of matrix A in Eq.(39) 



For simplicity, we put SH^H = a , = *. J = c, f = d, § = e, then 

A = i{{0, -e, -d, 0, e, 0, 0, 0, d, 0, 0, 0, 0, 0, 0, 0}, 

{-e, a-y-b-2c,0, -d, 0, e, 0, 0, 0, d, 0, 0, 0, 0, 0, 0}, 
{-d, 0,a-y + b-2c, -e, 0, 0, e, 0, 0, 0, d, 0, 0, 0, 0, 0), 
{0, -d, -e, 2a - Zy, 0, 0, 0, e, 0, 0, 0, d, 0, 0, 0, 0}, 
{e, 0,0,0, -a+y + b + 2c, -e, -d, 0, 0, 0, 0, 0, d, 0, 0, 0), 
{0, e, 0, 0, -e, 0, 0, -d, 0, 0, 0, 0, 0, d, 0, 0), 
{0, 0, e, 0, -d, 0, 2b, -e, 0, 0, 0, 0, 0, 0, d, 0], 
{0, 0, 0, e, 0, -d,-e,a-y + b + 2c, 0, 0, 0, 0, 0, 0, 0, d), 
{d, 0, 0, 0, 0,0,0,0,-a + y-b + 2c, -e, -d, 0, e, 0, 0, 0}, 
{0, d, 0, 0, 0, 0, 0, 0, -e, -2b, 0, -d, 0, e, 0, 0], 
{0, 0, d, 0, 0, 0, 0, 0, -d, 0, 0, -e, 0, 0, e, 0}, 
{0, 0, 0, d, 0, 0, 0, 0, 0, -d, -e,a-b + 2c-y, 0, 0, 0, e], 
{0, 0, 0, 0, d, 0, 0, 0, e, 0, 0, 0, -2a + Zy, -e, -d, 0\, 
{0, 0, 0, 0, 0, d, 0, 0, 0, e, 0, 0, -e, -a - b - 2c + y, 0, -d), 
{0, 0, 0, 0, 0, 0, d, 0, 0, 0, e, 0, -d, 0,-a + b-2c + y, -e], 
{0, 0, 0, 0, 0, 0, 0, d, 0, 0, 0, e, 0, -d, -e, 0}} 
-DiagonalMatri X [{± ±, ± ±, ± ±, ± ± ± ± ±, ± ±, ± ±, ±}]. 
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